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Abstract 


A  method  has  been  suggested  to  compare  and  evaluate  vibration  environments  or  sinusoidal 
and  random  specifications  through  the  use  of  a  single  analysis  technique.  The  comparison 
methodology  utilized  in  this  analysis  could  be  significantly  implemented  and  exploited  in 
interpreting  spectrum  analysis  results  of  vibration  data  such  as  in  the  analysis  of  a  gas  turbine 
engine  of  a  modem  tank.  A  generalized,  multi-degree-of-ffeedom  (DOF)  lumped  parameter 
structural  system  model  was  used  to  allow  an  initial  evaluation  of  the  environments  such  as 
shock,  sine,  and  random  so  as  to  possibly  eliminate  a  detailed  separate  dynamic  analysis  for  each 
environment.  Results  indicate  a  possible  reduction  of  analysis  time  and  costs.  In  addition,  since 
both  sine  and  random  analyses  depend  heavily  on  modal  damping  assumptions  for  the  accuracy 
of  their  predictions,  the  simple  methodology  proposed  herein  should  prove  to  be  useful  and 
productive. 
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1.  Introduction 


Vibration  analysis  involving  designs  should  be  evaluated  for  a  variety  of  loads  in  order  to 
determine  the  severity  and  possible  failure.  Not  one  environment  (viz,  shock,  sinusoids,  or  random 
excitation)  allows  the  design  of  all  structural  components  and,  hence,  we  are  led  into  separate 
analyses  for  shock — ^random  or  periodic  motions. 

This  work  purports  an  approach,  although  quite  simple,  to  compare  and  evaluate  vibration 
environment  or  sinusoidal  and  random  specification  via  a  single  analysis  technique.  This 
comparison  methodology  will  prove  to  be  extremely  useful  in  interpreting  spectrum  analysis  results 
of  vibration  data  that  will  be  taken  of  the  gas  turbine  engine  of  a  typical  modem  tank,  especially 
pertaining  to  the  turbine  engine  diagnostics  (TED)  project  (Helfman,  Dumer,  and  Hanratty  1995). 

2.  Methodology 

In  general,  consider  an  q  (multi)  degree-of-ffeedom  (DOF)  lumped  parameter  structural  system 
with  mass  matrix  [MJ,  stiffness  matrix  [kg],  damping  matrix  [Qj],  and  column  matrix  of  external 
forces  If  (x.,  t)}.  Mj,  kg,  Qj,  and  Fj  are  expressed  in  the  w  coordinate  system.  Given  a  forcing 
function  in  the  form  as 

IF  (x.,  t)}  =  {P^  p  (x.)}f  (t),  (1) 

then,  the  differential  equations  of  motion  in  the  w  coordinate  system  take  the  matrix  form 

[m]{w}  +  [c]w  +  [k]w  =  {p  (x.)}f  (t).  (2) 

If  we  apply  a  coordinate  transformation. 
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(w)  =  [<})]  {r|}. 


(3) 


in  which  each  column  of  (f)  is  a  modal  column  of  the  system  and  {t|}  represents  the  normal 
coordinates,  then,  from  equation  (3), 

(w)  =  (4) 

and 

{w}  =  [<j)]{ti}.  (5) 

Substituting  equations  (3)-(5)  into  equation  (2)  and  premultiplying  by  the  transpose  of  [(|)], 
equation  (2)  becomes 

[♦f  *  [‘I’f  Mmiti)  [<l>r  M[<l>](>i)  =  PJ<l>f  tp(!Cj))f(t).  (6) 

From  the  orthogonality  relations  of  natural  modes,  it  follows  that 

[<l>f  [m]  m  =  [MJ.  (7) 

and 

[<t>fM[<l>]  =  =  [K_].  (8) 

Comparing  the  triple  matrix  product 

[<i>f  w  m  (9) 
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with  equations  (7)  and  (8),  this  product  will  result  in  a  diagonal  matrix  only  when  the  damping 
matrix  [c]  is  proportional  to  either  the  mass  matrix  [m]  or  the  stiffness  matrix  [k];  that  is, 


w  =  2P[M], 

(10) 

[c]  =  a[k]. 

(11) 

or 


=  (12) 

where  Cr  is  the  ratio  of  the  assumed  to  critical  damping  and  cOj  is  the  angular  frequency.  Using 
equation  (10)  in  equation  (9)  and  comparing  with  equation  (7), 

[<l>f  W  [<H  =  2P[MJ.  (13) 

The  substitution  of  equations  (7)  and  (8)  with  equation  (6)  results  in  a  set  of  n  decoupled  differential 
equations  of  motion. 

[MJlfi}  +  2P  [MJ  (f|)  +  [o>;]  [M^  In)  =  [((>]’■  ip  (Xj) }  f(t).  (14) 

The  rth  equation  of  equation  (14)  has  the  form 

M^n,  +  2PM, n,  +  “'M,n,  =  P.f<l>(r))  lp(x))f(t).  (15) 


Let 


{<l)(r)}  {p(x.)}  =  =  participation  factor.  (16) 

Dividing  through  by  in  equation  (15)  and  using  equation  (16)  results  in 
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(17) 


f)  +2Pfi  +  o)^r|  =  —  r  f(t). 

'r  'r  r  *r  r  ^  ^ 

r 

A  general  form  of  solution  to  equation  (17)  is  obtained  readily  by  the  use  of  Laplace  transform 
(Hurty  and  Rubenstein  1964;  Gardner  and  Barnes  1942).  For  the  case  of  shock  spectrum,  the 
solution  to  equation  (17)  becomes 


D  (t). 


in  which 


D  (t)  =  f  ‘  h(t  -  t)  f(T)  dt 

J  o 


=  f‘  — 

Jo  I 


=  sin  ^6)^^  (t  -  t)  f(T)  dt , 


where  D/t)  is  called  the  dynamic  load  factor. 

The  total  deflection  of  the  structure,  w  (x,t),  is  obtained  by  inserting  equation  (19)  into 
equation  (3).  Thus, 

n  r. 

w(x,t)=P/S  ^L_4,D.(t).  (20) 

i  =  1  0)2  M  .  ’  ‘ 

1  1 

Starting  again  from  equation  (17),  and  posing  the  problem  of  harmonic  or  sine  excitation. 


- 


f(t)  =  e 


Use  of  Laplace  transform  provides  the  steady-state  solution  given  by 
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^,(t) 


p  r 

o  r 


1  -  -5-  .i2c 

2  r 

(0^ 
r 


p  r 

H(0)f(t), 

oo^M 


Q 


0) 


f(t) 


(22) 


(23) 


where 


r  =  1,  2,  n;  H  (Q)  = 


1  - 


Q 

V  r/ 


+  i  2C 


Q 


G) 


(24) 


The  total  deflection  of  the  structure,  w  (xjt),  is  obtained  by  inserting  equation  (22)  into  equation  (3). 
Thus, 


n  r. 

W(x,t)  =  S  <|>|  H,  (Q)  f(t), 

i  =  1  (*)"  M. 

1  i 


(25) 


or 


w(x,t)  =  P^  S  (j>^Ti^(t). 

r  =  1 


(26) 


The  mean  square  of  the  response  can  be  written  as 


1  T 


w  (x,t)  =  lim  —  J  w^(x,t)  dt. 

X-OO  -X 


(27) 


Substituting  in  equation  (27)  from  equation  (26)  and  then  from  equation  (23),  we  have 


w 


(x,t)  =  lim  c{)^  (x)  (j)^  (x)  ~  (t)  (t)  dt 


p2  r  r 


^  ^  ;7  ^  H  (Q)  (t)  dt. 

r  =  1  s  =  1  0)^  <0^  M  M  T-~  2T  J-T 


r  s  r  s 


(28) 
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If,  as  an  approximation,  we  disregard  phase  relations  that  will  tend  to  result  in  a  higher  mean  square 
value  in  equation  (28)  (i.e.,  neglecting  terms  that  are  a  significant  departure  from  the  mean 
[anomalies]),  then  the  integral  on  the  right-hand  side  of  this  equation  can  be  written  as 

Um  4;  f  ’■  H  (Q)  H  (Q)  f  ^  (t)  dt.  (29) 

X-co  •'-T 

When  the  forcing  function  f(t)  is  a  representative  record  of  an  ergodic  process,  we  can  transform 
the  limiting  process  of  equation  (29)  from  the  time  domain  to  the  frequency  domain  because  the 
function  f(t)  is  then  represented  by  frequency  components  in  a  continuous  frequency  spectrum 
0<Q<oo.  Thus, 


f\t)  =  liin  i  =  ^  ff(Q)dQ. 

x-«.  2T  J  -T  271  o 


Using  this  transformation  from  the  time  to  the  frequency  domain  in  equation  (29),  we  write 


Um  ^  |h,  (q|h  (Q)|f^  W  *  =  ^  /■  K  Hh.  H*'  ® 


Substituting  from  equation  (28),  we  write  for  the  mean  square  response  of  an  n  DOF  system  excited 
by  a  random  forcing  function: 


w^(x,t)  -  I,  Z,  d>.  X  ^  /:  K  (4.  (!3)|f  (D)  dQ,  (32) 


r  -  1  e  *  1 


0)^  CO^  M  M  271  Jo 


The  solution  of  equation  (32)  can  be  simplified  further  assume  a  lightly  damped,  multi-DOF  system. 


H  (Q)  = 


QV  ,r2 

—  +  4C  — 

«  ^  0)2 
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According  to  Hurty  and  Rebinstein  (1964),  the  magnification  factors  (Q)  have  regions  of 
pronounced  peaks  in  the  neighborhood  of  the  corresponding  natural  frequencies  w,.  Furthermore, 
Biggs  (1964)  tells  us  that  the  products  (Q)  and  (Q)  for  r  #  s  are  seen  to  be  small  in 
comparison  with  the  same  products  for  r  =  s.  In  addition,  in  equation  (32),  terms  with  r  s  may  be 
negative  as  well  as  positive,  depending  on  the  sign  of  the  product  (j)^  (j)^  F^ ,  while  terms  with 
r  =  s  are  always  positive.  Equation  (32)  will  therefore  be  small.  Hence,  by  disregarding  the  cross- 
product  terms  corresponding  to  r  =  s,  equation  (32)  becomes 

w"(x,t)  =  I  <j>^^  i  r  |h  (Q)  (Q)  dQ,  (34) 

r  =  1  (0^  2lt  Jo  I  " 

r  r 

Equation  (34)  can  be  further  simplified  by  approximating  the  integrals  of  equation  (34)  by  replacing 
f  (Q)  by  its  discrete  values  f((0r)  at  the  natural  frequencies  Wj  and  using 


and 


(36a) 

S'i.  = 

(36b) 

f  =  f  =  27ttO  , 

nr  r’ 

(37) 

(38) 
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equation  (34)  becomes 


w  (x,t)  =  Z  %  (^) 

r  =  1 


p2r2  ( 

o  r 


Tl  ,  _  _ 

—  g  .  f  Q 

2  ®  m  n  ^ 


(39) 


Equation  (39)  is  the  equation  that  forms  the  basis  of  the  method  of  response  spectrum  analysis. 
Hence,  in  this  perspective,  we  have  restricted  our  analysis  to  the  following  constraints: 


•  Ignored  phase  relations. 


•  Ignored  magnification  products  at  different  frequencies. 


•  Lightly  damped  all  modes  (frequencies),  and 

•  Used  discrete  values  of  the  input  spectrum  at  the  natural  frequencies. 

These  conditions  are  realistic  if  we  are  dealing  with  a  system  that  is  grounded  by  base  mounts,  and 
that  has  harmonic  frequencies.  At  this  point,  we  can  summarize  and  rewrite  equations  (20),  (25), 
and  (39)  as  follows: 


•  Case  no.  1 — shock  (from  equation  [20]): 


w(x,t)  = 


P 


O 


r 

r 

0)2  M 

r  r 


«{>,  (x)  (t), 


where  r  =  i  =  1,  2, ...,  n;  or  from  equation  (25), 


n  p  r 

W(x,t)  =  S  (|)^  (x)  d;  (t). 


1  co^  M 


where 
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dt^ 


*  Case  no.  2 — sine  (from  equation  [39]): 

w  (x,t)  =  J  -  I  - "  (j)^  (x)  H  (Q)  f(t), 

r  =  1  0)2  " 


where  r  =  i  =  1,  2, ...,  n;  or 


w(x,t) 


n  p  r 

S  ({>^(x)Qf' (t), 

r  =  1  0)2  M 

r  r 


where 

f  (t)  =  -^  («t)), 

dt^ 

and 


(from  equation  [33]),  when  Q  =  0),. 


•  Case  no.  3 — ^random: 


w^(x,t)  =  X 


.  p2r2 

o  r 

1 


/ 


(x) 


7Z  .  _ 

—  g  f  Q 

2  ®  m  n  ^ 


(40) 


where,  in  all  cases,  w(x,t)  is  defined  as  response  acceleration.  In  case  no.  3,  response  gn^s 
is  defined  as 


Srms 


=  Vw"(x,t)  = 


p2r2 

o  r 

(O*  m2 

r  r 


4>J  (x) 


'f  °  m  n  ^ 
^  ) 


1/2 


(41) 
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which  includes  all  frequencies  (w„  Q  and  modes  (|>Xx)  of  the  structural  system  under 
consideration.  If  we  look  at  gR^g  per  mode  or  frequency,  we  are  led  to 


p  r 

o  r 

0)2  M 


r  « 


TZ 


f  Q 


1/2 


(42) 


as  a  first  step  toward  the  overall  gRMs  computation.  By  comparing  equation  (41)  to  case  no.  2  and, 
in  turn,  to  case  no.  1,  we  see  that  the  only  difference  in  the  response  of  a  multi-DOF  structure  at  a 
natural  frequency,  Wj,  are  the  forcing  function  terms;  that  is, 

case  no.  1 — shock; 


d;  (t). 


(43) 


case  no.  2 — sine: 


Q  f'  (t). 


(44) 


and  case  no.  3 — ^random: 


/ 


7t 

2 


g'  f  Q 

®  m  r  ^ 


in 


All  other  terms  of  the  aforementioned  cases;  that  is. 


(45) 


p  r 

~~  (x),  (46) 

M 

r  r 

define  a  multi-DOF  system  (as  opposed  to  a  one-DOF  system)  and  are  common  to  all  gresponse 
solutions  for  shock,  sine,  or  random  environments. 
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Equations  (42)-{44)  are  nothing  more  than  the  gresponse  terms  as  discussed  for  a  simple  one-DOF 
system  by  Hurty  and  Rubinstein  (1964). 

Using  the  analysis  technique  for  response  spectrum  analysis  and  saving  the  gn„s  terms  at  each 
natural  frequency  will  then  allow  the  determination  of  true  gnns  response  for  each  environment 
(shock,  sine,  or  random)  that  differs  by  a  gi„pu,  magnitude  difference  defined  and  known  via 
equations  (42)-(44). 

3.  Conclusions 


The  approach  shown  previously  provides  the  benefits  of  time  and  cost  savings  but  is  limited  by 
the  assumption  required  to  derive  the  approach.  The  application  of  the  approach  can  be  verified 
early  in  the  design  but  requires  the  capability  to  solve  the  dynamic  analysis  using  a  program  that 
utilizes  both  approaches  formulated  by  equation  (28)  vs.  equation  (39). 

Notwithstanding  the  fact  that  shock  response  is  little  affected  by  assumed  damping,  both  sine 
and  random  analyses  depend  heavily  on  modal  damping  assumptions  for  the  accuracy  of  their 
predictions.  Although  the  approach  in  this  analysis  was  based  upon  the  concept  of  response 
spectrum  analysis,  its  concept  is  used  quite  extensively  in  oscillatory  systems  and  even  in  earthquake 
analysis.  Here,  the  simplified  equations  can  be  directly  correlated  to  sine  or  shock  vibration 
analyses.  Hence,  by  means  of  this  simplified,  but  yet  practicable  technique,  one  can  easily  solve  all 
or  any  one  of  the  designing  environments  for  a  given  structure,  such  as  the  case  for  the  mechanical 
structural  components/mechanisms  inherent  in  a  typical  gas  turbine  engine  assembly  network. 

However,  this  approach  is  feasible  provided  that  phase  relations  are  ignored  and  magnification 
factors  are  not  a  problem.  Furthermore,  this  methodology  is  especially  useful  for  systems  involving 
harmonic  frequencies  based  upon  the  natural  frequencies  of  the  individual  mechanical  components. 
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